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Abstract. We describe the transfer of quantum information and entanglement from 
three flying (radiation) to three localized (atomic) qubits via cavity modes resonantly 
coupled to the atoms, in the presence of a common reservoir. Upon addressing the 
full dynamics of the resulting nine-qubit open system, we find that once the cavities 
are fed, fidelity and transferred entanglement are optimal, while their peak values 
exponentially decrease due to dissipative processes. The external radiation is then 
turned off and quantum correlations oscillate between atomic and cavity qubits. For 
a class of mixtures of W and GHZ input states we deal with a discontinuous exchange 
of entanglement among the subsystems, facing the still open problem of entanglement 
sudden death and birth in a multipartite system. 
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1. Introduction 



As early as in 1935 Einstein with Podolski and Rosen [TJ and Schrodinger [2] drew the 
attention on the correlations in quantum composite systems and the problems raised by 
their properties. Much later, theoretical [3] and experimental [4j cornerstones elucidated 
the issue of nonlocality. Entanglement is currently viewed as the key resource for 
quantum information processing [5], where it allowed a number of achievements such as 
teleportation [6], cryptography [7] and enhanced measurements [H]. The deep meaning 
of multipartite entanglement, its quantification and detection [9] as well as its possible 
applications, are the object of massive investigation touching the heart of quantum 
physics and a great variety of its branches. In this paper we study the entanglement 
dynamics of a multipartite open system in cavity quantum electrodynamics (CQED) 
[TO] by a model that is not so far from physical implementation with current technology 
of optical cavities [UJ. In particular we study a paradigmatic example to investigate 
multipartite entanglement transfer and swapping [T2J, IT5] . which are fundamental 
processes for the implementation of quantum interfaces and memories for quantum 
networks [HI EE] • In addition we deal with the recently discovered [16] and observed [T7] 
phenomenon of entanglement sudden death (ESD) (and birth (ESB)), consisting in an 
abrupt vanishing (and raising) of quantum correlations. The problem is nowadays still 
open with respect to the multipartite entanglement classification, but we can shed some 
light at least on the ESD and ESB of the fully tripartite entanglement. The interplay 
among all these aspects in the presence of external environments is investigated. 



2. The multipartite open system model 



We consider three entangled radiation modes which are injected and resonantly cou- 
pled with three separated optical cavities, each of them containing a trapped two- 
level atom whose transition frequency is resonant with the cavity mode. Each cav- 
ity mode is coupled to an external environment by an amplitude damping channel 
with a decay rate k and each atom can spontaneously emit a photon with a decay 
7. The evolution of the whole system density operator p(t) can be described by a 
Master Equation (ME) in the Lindblad form that we have solved by means of the 
Monte Carlo Wave Function method [18J. We identify a set of six collapse operators 
C c ,j = Vhcj and C a ,j = \Z^o~j (J=A,B,C), where k = k/g a and 7 = 7 / 'g a are di- 
mensionless cavity and atomic decay rates scaled to the atom-cavity coupling constant 
g a (taken equal for each cavity mode-atom subsystem). The effective Hamiltonian is 
Tie = y~ — y J2j &l jC c j + C\ jC a: j , with the system Hamiltonian in the interaction 

picture H 1 = J2j [^9a(cja\ + c\aj) + ihg c (cjf\ - c ] jfj)}. Here cj, c]j (fj, f\) are the 
quantum harmonic oscillator operators for the cavity (external radiation) modes, &j, a\ 
the raising and lowering operators for the atomic qubits, and g c the cavity-input field 
coupling constant. 
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2.1. Subsystems dynamics, state mapping and entanglement transfer 

We have solved the dynamics for several types of entangled external radiation prepared 
in mixed or pure states, cavity mode prepared in the vacuum state and different types of 
atomic preparation. Here we concentrate on the case of atoms in the lower energy state 
\ggg) a and the external field prepared in a qubit-like entangled state because this is the 
condition for maximum entanglement transfer to the atomic subsystem [13]. Overall 
we are thus dealing with a 9-qubit system. From now on, for sake of simplicity, we 
take equal coupling constants of the cavity modes with the atoms and the external 
radiation (g c = g a = g) and we deal with dimensionless times r = gt. We first consider 
the case of negligible cavity and atomic decays. We show that it is possible to map 
the initial external radiation states onto the atomic and the cavity subsystems states, 
thereby also transferring the initial quantum correlations. Due to the qubit-like form of 
the external and cavity mode fields we quantify the tripartite subsystems entanglement 
by the tripartite negativity E^(r)(a = a,c,f) [19] . defined as the geometric mean 
of the three bipartite negativities [20] . In Fig. [T] we illustrate results for the case of 
external radiation prepared in the W state |tf(0))/ = (|001)/ + |010)/ + |100> / )/v / 3- 
We notice that each subsystem (A,B,C) has the same dynamics that is composed by 
a transient and an oscillating regime. In the transient each flying qubit transfers its 
excitation to the cavity qubit which in turn passes it onto the atomic one, as shown in 
Fig. QJi. This is the physical process for the transfer of quantum information. The cavity 
mode, simultaneously coupled to the external field and the atom (describable as a third 
harmonic oscillator), exchanges energy according to a Tavis-Cummings dynamics at an 
effective angular frequency gy2 [TDj, [21]. During the transient up to time r ff = ir/y/2 
the mean photon number N^(t) = ($c)(t) in each cavity completes a cycle. In 
the same period the atomic excitation probability p e (r) reaches its maximum value, 
while the input field has completely entered the cavity, i.e. its mean photon number 
NW(t) = (Pf)(r) vanishes. 




Figure 1. Dynamics for the external field in a W state: (a) mean photon numbers 
(dashed), N^' (dotted) and atomic probability p e (solid); (b) tripartite negativity 
ij( Q ) f or atoms (solid), cavity modes (dashed) and external field (dotted); ( c) atomic 
fidelity (solid) and purity /^ a ) (dashed); (d) cavity mode fidelity F( c ) (solid) and 
purity fi^ c \ 

In Fig. [lb we show and see that up to r ff/2 the entanglement of the external 
driving field is mainly transferred to the cavity modes and then allows the build up of 
atomic entanglement. We note that E^(r ff) = 0.94 which is the value of the tripartite 



Tripartite entanglement transfer in a CQED open system 



4 



negativity for the injected W state. For the time evolution of the atomic and cavity mode 
subsystems described by the density operator p a (r) (a = a,c) we show in Figs. H](c,d) 
the purity = Tr Q [p2( r )] and the fidelity F^{t) = a (*(0)|p a (r)]*(0)) Q with 

respect to a pure state |\1/(0)) Q that is a map of the initial external field 1^(0))/, with 
the correspondence |0)/ <-» \g) a (|0) c ) and <-» \e) a (|l) c ) for the atoms (cavity 
mode) states. We see that up to r ff the dynamics maps the whole initial state 
|*(0))/®|000) c ® \ggg) a onto |000)/® |000) c ® E/_»|tf(0)) a , where fjf = is 
a local phase rotation, that in the case of W states simply acts as a global phase factor 
+i. Note that the maximum of E^ c '(r ff/2) does not correspond to a pure state, i.e. 
the initial state |\l/(0))/ cannot be mapped onto the cavity modes during the transient 
regime. 

At the end of the transient regime the external radiation is turned off and the 
subsequent dynamics can be described by a triple Jaynes-Cummings [22] ruled by 
oscillations at the vacuum Rabi frequency 2g, hence with a dimensionless period tt 
as shown by cavity mean photon number and atomic probability in Fig. [Th- The 
purities /^ a > c ) (r) in Figs.[T](c,d) oscillate at a double frequency between pure full tripartite 
entangled and separable states. In particular, at times r m = r ff + mir (m = 0, 1,2...) 
the atoms are in the entangled states C/^|^/(0)) a , where = =p| and — (+) applies in 
correspondence to even (odd) values of m, that are the peaks of £™(t) in Fig. [lb. At 
times r n = r ff + (n + |)7r (n = 0, 1,2...) the cavity modes are in the qubit-like state 
£7^ |\&(0)} c , where = e -i<i>cjcj j g a \ 0C8k \ rotation such that 4> = (<j) = re) for even 
(odd) values of n. As regards the bipartite subsystem entanglement, evaluated after 
tracing over all other qubits, we find sizeable two-atom and two-cavity entanglement 
(negativity peaks ~ 0.41), as well as some entanglement for all other qubit pairs (equal 
negativity peaks ~ 0.08). 

Starting, instead, from an external field in the qubit-like GHZ state |\&(0))/ = 
4|(|000)j + |111)/) we find a quite similar time evolution but different values for the 
maxima of the quantities evaluated in Figs. [T^a,b). In particular, at times r m we have 
the maximum entanglement transfer with = 1, that is again the same value of 
|^(0))/, and the atomic states are mapped onto (|000} o =F i|lll) a )/\/2. On the other 
hand, at times r n the cavity modes are in the state (|000) c ± |lll) c )/\/2. The two-qubit 
entanglement is null for any pair, reflecting the lack of subsystem entanglement in a 
GHZ state, except the three directly coupled atom-cavity pairs, where the negativity 
peaks occur at about 0.21 at any quarter and three quarters of a period. We remark that 
the above state mapping can be obtained also for any qubit-like state |\l/(0))/ written 
in a generalized Schmidt decomposition [23] . 

2.2. Dissipative effects 

In the perspective of experimental implementation of our scheme for quantum 
information purposes an important issue is the effect of dissipation on both state 
mapping and entanglement transfer. For cavity decay rates in the range < k < 0.4 
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we calculated the fidelities F^ a \r) at the first peaks (after the transient) as well as 
the tripartite negativities at the first two peaks E^ a ' of both atomic and cavity field 
subsystems, for external field preparations in states W and GHZ. The behavior of all 
these quantities can be well described as f{k) = /(0)e _/3fc , where the values of decay 
rates (3 are reported in Table [U As expected, quantum state mapping and entanglement 

Table 1. Effect of cavity decay rate k on F^ a >{T m ) at the first peak (m = 0) and 
E^(r m ) at the first two peaks (m = 0,1), and on F^(r n ) (n = 0) and £ (c) (r„) 
(n = 0, 1) for injected field states W and GHZ. 





F^(t ) 




£( a )(ri) 


F^(r ) 


E^(t ) 


£ (c) (n) 




0.55 


1.09 


4.59 


1.41 


3.22 


6.47 


Pghz 


0.77 


1.10 


4.69 


1.88 


2.80 


6.11 



transfer are by far more efficient onto atomic than cavity qubits. Actually these processes 
involve the atoms when the cavity modes, coupled to the external environment, are 
negligibly excited. Cavity modes are instead affected by dissipation already in the 
transient regime, which leads to a reduction of the relevant peaks as well as asymmetries 
in their profile (decay slower than build up). All these features, imprinted in the 
transient, become quite manifest in the subsequent cavity dynamics. 

We also analyzed the effect of the spontaneous emission of atomic excited levels 
on state mapping and entanglement transfer processes. We consider, as an example, 
the case of decay rate k = 0.1 and values of dimensionless parameter 7 up to 0.1. The 
exponential approximation is less accurate and we report in Table [2] the results for 
exponential decay rate 5 only at the first peaks of fidelity and tripartite negativity for 
both atomic and cavity mode subsystems. 

Table 2. Effect of atomic decay rate 7 on F^(T m ) and E^ a >{r m ) at the first peak 
(m = 0), and on F^(T n ) and E^(r n ) for n = 0, for injected field states W and GHZ. 





F^(t ) 


£ (a) (To) 


FW(ro) 






1.01 


2.06 


1.85 


4.07 


Sghz 


1.57 


2.21 


3.14 


4.51 



3. Fully tripartite ESD and ESB 

Let us now consider the class of mixed qubit-like states for the injected field p/(0) = 
p\GHZ)(GHZ\ + (1 — p)|W / )(W|, < p < 1. By analytical and numerical results we 
show that our scheme for entanglement transfer and swapping is also relevant for the 
observation of sudden disentanglement and entanglement effects. Since the tripartite 
negativity E^\t) with (a = a,c) is an entanglement measure that provides only a 
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sufficient condition for entanglement detection, we cannot properly talk about ESD 
or ESB for these kind of states in the whole parameter space {p, r}. Nevertheless, 
we can classify the atomic state by using the entanglement witnesses |24J Wg = 
| J - \GHZ)(GHZ\, W W 2 = \I ~ \GHZ}(GHZ\ and W W i = %I - \W)(W\, and 
analyze the discontinuous evolution of entanglement focusing only on the fully tripartite 
entanglement properties. For negligible dissipation and at times r m , the initial state of 
radiation is mapped onto the atoms, then the entanglement classification is known. 
The state is of class GHZ for | < p < 1, of class W for < p < | and 
\ < V < 1 > an d biseparable (B) for | < p < ~. Outside times r m , we can still 
make a partial entanglement classification that is shown in Figj2^L (Figj2b) for atomic 
(cavity field) states. Since the tripartite negativity is zero in the black regions of FigfSl 
we cannot exclude the presence of bound entangled states. With this knowledge on 
the entanglement properties in parameter space, we can anyway affirm that there is 
a discontinuity for the full tripartite entanglement. Fixing, for instance, a value of 
p < 0.25 and looking at the time evolution of the atomic state, we can notice that it 
suddenly acquires a fully tripartite entanglement entering the W region and loses this 
property after some finite time exiting that region. This effect can be addressed as an 
ESD and ESB of the fully tripartite inseparability only. Moreover, since the systems 
share among its subsystems energy and entanglement in a periodic way, we can highlight 
also that discontinuities in the fully tripartite entanglement are exchanged among them. 
This occurs after the transient, where the cavity modes do not exhibit genuine tripartite 
entanglement (see Fig. [2b) but only support the transfer of quantum correlations from 
the input field to the atoms. 




Figure 2. Entanglement classification in the parameter space {p, r} for external field 
prepared in the mixed GHZ and W state /5/(0). a) atomic subsystem, b) cavity mode 
subsystem. 

4. Conclusions 

In this paper we have addressed the transfer of quantum information and entanglement 
from three flying qubits to three localized ones, focusing on the basic physical features 
characterizing multiqubit state mapping in a CQED setting. We analyzed the effect of 
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dissipation at the times when the transfer protocol is optimal, considering the atoms 
and the cavities both in contact with a same reservoir at zero temperature. We derived 
also the conditions for the repeated occurrence of discontinuous exchange of quantum 
correlations among the multipartite subsystems. Our scheme could be implemented by 
combining current advances with optical fibers, cavities, and trapped atoms [TT1 125| 126]. 
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